Physical Realizability of Multi-Level Quantum Systems* 
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Abstract — This paper considers the physical realizability 
condition for multi-level quantum systems having polynomial 
Hamiltonian and multiplicative coupling with respect to several 
interacting boson fields. Specifically, it generalizes a recent 
result the authors developed for two-level quantum systems. 
For this purpose, the algebra of SU(n) was incorporated. As 
a consequence, the obtained condition is given in terms of the 
structure constants of SU(n). 

I. Introduction 

In an environment where the classical laws of physics 
apply, standard control techniques such as optimization or a 
Lyapunov procedures do not worry in general of the nature of 
the controller they synthesized. In other words, their imple- 
mentation is always possible since the physics behind them 
still hold. However, if one desires to implement a controller 
that obeys the laws imposed by quantum mechanics (e.g., 
quantum coherent control [2], [7], [12]), then such a task 
is not so easily achieved unless an explicit characterization 
of those laws is given in terms of the control system vector 
fields. This is exactly the purpose for introducing the concept 
of physical realizability. 

Conditions for physical realizability were first given 
specifically for linear systems satisfying the quantum har- 
monic oscillator canonical commutation relations [5], [8]. 
Recently, the formalism was extended for systems describing 
the dynamics of open two-level quantum systems interacting 
only with one quantum field in which the algebra of SU(2) 
played a central role [3]. Compared to a linear quantum 
system, the systems being analyzed were bilinear, and the 
the commutation relations were dependent on the system 
variables. Thus, the main contribution of this paper, given in 
Section [IV] is to provide a condition for physical realizability 
of multi-level quantum systems having polynomial Hamilto- 
nian and multiplicative coupling, and whose system variables 
obey the commutation relations described by the algebra 
of SU(n). As expected, the obtained condition is given in 
terms of the symmetric and antisymmetric structure constants 
of SU(n). Another contribution is that the systems under 
consideration have been allowed to interact with multiple 
quantum fields in quadrature form. 

The paper is organized as follows. Section HH presents the 
basic preliminaries on open quantum systems. In Section HIT] 
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the necessary algebraic machinery to study open multi-level 
quantum systems is given. This is followed by Section [IV] 
in which the definition of physical realizability is provided 
as well as a condition for a bilinear QSDE to be physically 
realizable. Finally, Section IVl gives the conclusions. 

II. Open Multi-Level Quantum Systems 

Open quantum systems are systems governed by the laws 
of quantum mechanics that interact with an external environ- 
ment. A quantum mechanical system is described in terms 
of observables and states. Observables represent physical 
quantities that can be measured, as self-adjoint operators 
on a complex separable Hilbert space Sj, while states give 
the current status of the system, as elements of f), allowing 
the computation of expected values of observables. In [1], 
[11], the evolution of open quantum systems is given in 
terms of quantum stochastic differential equations. For this 
purpose, observables may be thought as quantum random 
variables that do not in general commute. A measure of the 
non commutativity between observables is usually given by 
the commutator between operators. 

Definition 1: The commutator of two scalar operators x 
and y in fj is the antisymmetric bilinear operation 

[x, y] = xy- yx. 

Also, if a; is an m-dimensional vector of operators in and 
y is an rii -dimensional vector of operators in fj, then 

„T\T 



r Ti A T 

[x,y \= xy 



(yx 



which is an m x n% matrix of operators in ¥). 
This commutator satisfies 



,TiT 



l%,y* r =-v x + ( x y ) =-[y,x } 

The adjoint of x is denoted by x' = (x#) T with 



(1) 



V < J 

and * denotes the operator adjoint. In the case of complex 
vectors and matrices, * denotes the complex conjugate while 
t denotes the conjugate transpose. The non-commutativity 
of observables is a fundamental difference between quan- 
tum systems and classical systems in which the former 
must satisfy certain commutation relations originating from 
Heisenberg uncertainty principle. 

The environment consists of a collection of oscillator sys- 
tems each with annihilation field operator w(t) and creation 
field operator w*(t) used for the annihilation and creation 
of quanta at point t, and commonly known as the boson 



quantum field (with parameter t). Here it is assumed that t 
is a real time parameter. The field operators w(t) and w*(t) 
satisfy commutation relations as well. That is, 

[w(t),w*(t')] =S(t-t') 

for all t, t' G K, where S(t) denotes the Dirac delta. Its 
mathematical description is given in terms of a Hilbert space 
called a Fock space. When the boson quantum field is in 
the vacuum state, i.e., no physical particles are present, it 
then represents a natural quantum extension of white noise, 
and may be described using the quantum Ito calculus [1], 
[11]. This amounts to having three interacting signals (inputs) 
in the evolution of the system: the annihilation processes 
W(t), the creation process W*(t), and the counting process 
A w (t). 

In simple words, the evolution of an open quantum system 
is described putting together the evolutions of the system and 
the environment in an unitary fashion. That is, if ip is an 
initial state then ip(t) = U(t)tp, where U(t) is unitary for 
all t, and is the solution of 

dU(t) = ^(S-I) dA w (t) + L dW* (t) -L*S dW(t) 

1 



{L*L+iH)dtj U(t), 

with initial condition U(0) = I, I denoting the identity 
operator and i being the imaginary unit. Here, % is a fixed 
self-adjoint operator representing the Hamiltonian of the 
system, and L and S are operators determining the coupling 
of the system to the field, with S unitary. The evolution of 
ip is equivalent to the evolution of the observable X given 
by 

X(t) = U*(t)(X®I)U(t), 

whose evolution is referred as the Heisenberg picture while 
the one for ip is known as the Schrddinger picture. This 
paper exclusively takes the point of view of the Heisenberg 
picture. The quantum stochastic calculus in [11] allows to 
express the Heisenberg picture evolution of an scalar operator 
X interacting with a boson field as 

dX = (S*XS - X) dA w + C(X) dt + S* [X, L] dW* 
+ [L*,X]SdW, 
where C{X) is the Lindblad operator defined as 

C(X) = -i[X, U] + \ (L* [X, L) + [L\X]L) . 
The output field is given by 

Y(t) = U(t)*W(t)U(t), 

which amounts to 

dY = Ldt + SdW. 

In summary, the dynamics of an open quantum system is 
uniquely determined by the triple (S,L,H). Hereafter, the 
operator S is assumed to be the identity operator (S = I). If 
on the other hand one consider n w interacting boson fields 
then the evolution equation is written as 

dX = C{X) dt + dW^ [X, L] + [L\X] dW, 



where [X, dW] = [X, dW^] T = 0, L = (L u . . . , L S ) T , 
C(X) ^ -i[X, H]+ l - (it [X , L] + [L\X] L) , 

dWx \ 

: and dW^ = (dW* , • • ■ , dW* w ) 



Consider the vector of operators x = (xi, . . . ,x s ) T . By 
stacking (column-wise) the scalar evolutions for each Xi, it 
follows that the Heisenberg evolution equation is 




dW* 




dW* 



= C(x) dt + [x, L T ] dW* - [x, L f ] dW, (2) 

where 

C(x) 
= -i[x,H] 

( tf[xi,L] \ ( \L\x x \L 



\ P,x 8 ]L 



i[x,H\-^(jL^([L, Xl ],-'- ,[L,x t ])) T 
\[x,tf]L 

i[x,H] + \(-(tf [L,x T }f -[x,L*]L) 

i[x,H] + \ ((^ [x,L T f) T + [L#,x T fL) . (3) 



It is customary to express QSDEs in terms of its interaction 
with quadrature fields. The quadrature fields are given by the 
transformation 

Wl\ = ( In w In u , \{W 

w 2 J \ -il nw il nw )\w* 

where the operators Wi and W2 are now self-adjoint, and 
I nm denotes the identity matrix of dimension n w . In [4], the 
Ito table for W and is 

Z* ) ( ^ ) - ( I '7 



(4) 



dt 



which in terms of the quadrature fields is 
Thus, 

dx = C(x) dt 



-iT T 



dt. 



1 



m ( I nw -il nw \ ( dW 1 
dW 2 



dt 



dWi 
dW 2 



([ Xi l t ],-[ x ,l^)( y 

= C(x) dt+- ([x, L T ] - [x, dW x 

- % -(\x,L T ] + [x,L^])dW 2 

The quadrature form of the output fields is obtained from 
the quadrature transformation 

Yl \ ( In w In w \(Y 
% J \ -Un w Un w A yt 

which gives 

d% \ f L + L* 
d% J { i{L# - L) 

The main focus of this paper is on the dynamics of 
open multi-level quantum systems interacting with n w Boson 
quantum fields. Such systems evolve with respect to the 
group SU(n). The algebra of SU(n) has been extensively 
studied since the 1950's to the point that it is an standard 
topic in quantum mechanics when studying multi-level sys- 
tems [9], [10]. To particularize the framework presented in 
the previous paragraph for system (f2]i evolving on SU(n), 
consider the Hilbert space S) = C" and let |j) with j = 
1, . . . , n be eigenvectors spanning S). A projection operator 
Pki is defined as the outer product 

Pk,i = \k) 

where k,l ~ 1, ...,n. It is a well-known fact that any 
operator defined in S) can be obtained in term of these n 2 
projection operators. Specifically the generators of SU (n) 
are constructed as follows 

u jk = Pj,k + Pk,j, 

v jk = i (Pj, k - P k!j ) , 



Wl 



1(1 + 1) 



P s ,s - kP i+u+ i 



for 1 < j < k < n, 1 < I < n— 1. The (n 2 — n)/2 symmetric 
matrices Ujk, the (n 2 — n)/2 antisymmetric matrices Vjk and 
the n — 1 mutually commutative matrices wi together form 
the set {Ai, . . . , A,^^]^} of generators of SU(n). Hereafter 
define s = n 2 — 1. Their commutation and anticommutation 
relations are 



[A;, Xj] — 2i ^2 fijkXi 



k=l 



{A;, Xj} = -Sij + 2 d ijk X k - 

k=l 



Thus, the product Ai Aj can be easily computed as 

XiX^^dX^XA + iXi,^}) 



2 

= —Sij + (ifijk + dijk) Xk 



(5) 



k=l 



where the real completely antisymmetric tensor /«/. and the 
real completely symmetric tensor duf. are called the structure 
constants of SU(n), and Sij is the Kronecker delta. The 
tensors fak and dijk satisfy 

filmfmjk ~t~ fjlmfimk ~t~ fklmfijm 0, (6) 



filmdmjk ~t~ fjlmdimk ~t~ fklradijm — 0: 



and 



^ . fimkfi 



jrnk 



nSi 



(7) 



(8) 



,k=l 



The procedure of how to construct the generalized Gell- 
Mann matrices shows that only the dimension of the group 
SU(n) is necessary to express all the components of the 
group algebra, i.e., once n is given then the generators and 
structure constants / and d are fixed. The vector of system 
variables for (f2j) is 

x\ \ ( Ai 

where Ai,...,A s are spin operators. Given that these op- 
erators are self-adjoint, the vector of operators x satisfies 
x = x# . In particular, a self-adjoint operator A in S) is 
spanned by the generalized Gell-Mann matrices [10], i.e., 

1 



A 



-a 



1 s 

It, 



oaXi, 



i=0 



where ao 



Tr(A), a t 



Tr(AAi). Thus, ao and 



(ai, . . . , a s ) G C determine uniquely the operator A with 
respect to a given basis in C™. The initial value of the system 
variables can be set to a;(0) = (Ai, . . . , A s ). 

It is important to emphasize that any arbitrary polynomial 
of the components of x evolving in SU (n) is described by 
a linear combination of the generators of SU(n), which 
includes the identity [10]. This fact can be appreciated in 
relation (0 because equation (f2]i implies that such com- 
mutation relations are preserved by the evolution in time. 
Thus, any Hamiltonian and coupling operators of polynomial 
type are representable as linear functions of x. Therefore, 
assuming linearity captures a large class of Hamiltonian and 
coupling operators without much loss of generality, i.e., the 
assumed Hamiltonian is TL — ax with a G M. s , and the 
multiplicative coupling operator is of the form L = Ax 
with A G iC n ™ xs . The reason why a coupling operator is 
called multiplicative is that they make the interacting fields 
to appear in (O as multiplicative quantum noise. 

In general, the evolution of x in quadrature form falls into 
a class of bilinear QSDEs expressed as 

dx = Agdt + Ax dt 

+ (BxiX, ■ ■ ■ ,Bx nm X, B 2 lX, ■ ■ ■ ,P>2n w x) ( J, (9) 



(10) 



where Aq G R s and A, B lk = B lfe + B 2fc , B 2fc 4 i(£ 2fc - 
Bik) G R sxs , fc = 1, . . . ,nu). The fact that all matrices in 
(0 are real is due to the quadrature transformation The 
quadrature output fields are whose quadrature form is 

dY 1 \_[C 1 \ ( dWi 

d% J ~ \ c 2 ) xm +y d w 2 

where C G C"«" xs , C\ = C + C* and C 2 = i(C# - C). 
Note that the quadrature transformation makes C\ and C 2 to 
be real matrices. 

The objective of this paper is to determine conditions on 
the coefficients in (0 and ( TTOb under which there exists the 
Hamiltonian and the coupling operator of the above form 
such that (0, ( [Tol l can be written as |2j see Definition [3] in 
Section |IV] 

III. Notation and Preliminary Results 

Define F,A e R sxs , i G {l,...,s}, such that their 
(j,k) component is {Fi) jk = f ljk and (A)jfc = <Ujk, 
respectively. In particular, the set {— iFx, . . . , — iF s } is the 
adjoint representation of SU(n). In [6], [9], identities (O 
and (0 were employed to obtain the following useful rela- 
tionships 



F i D j 



DiFj 



fijkFk 

[Fi,Dj] = -J2fijkD k 
k 

s 

/] dijkFk 



FjDi 



(11) 



(12) 



(13) 



DjFi 



k 

s 



— > dijkFk. 



k 



Definition 2: Let j3 G C s 
s — ► C sxs are defined as 



(14) 

The linear mappings Q~ , Q + : 



9-(/3) = (^ T /3,--- ,Fj/3) 



9+(P) = (D*0,~-,DT0) 



V P T FI 



(15) 



(16) 



In order to simplify the notation, if (3 is an s-dimensional 
row then it will be understood hereafter that _ (/3) = 
e~(/3 T ) and 9+(0) = 0+(/3 T ). In addition, these two 
matrix functions are used to express the commutation and 
anticommutation relations of the vector of operators x in a 
compact form. That is, 

[x,x T ] = 2iQ-(x), (17) 

{x,x T } = -7 + 26 + (x). (18) 
n 

It is also important to notice that the nature of the / 
and d-tensors make the matrices 0~(/3) and + (/3) be 
antisymmetric and symmetric, respectively. Consider now the 
stacking operator vcc : C mxn — > C" m whose action on a 



matrix creates a column vector by stacking its columns below 
one another. With the help of vec, the matrices 0~(/3) and 
+ (/3) can be reorganized so that 

/ er(/?) 

vec(6-(/3)) = : 



and 



/ e+(/3) 



vcc(0+(/3)) 



where /3 G C s , 0r(/?) = if/3, F = (Fi, • • • , F s f, 
Q+(/3) = D t j3 and D = {D u --- ,D s f. From ®, F 
satisfies 



Vfc=l / ij 

s 

^ {Fk)im{Fk)r 



k,m—l 



^ ^ fkim fkmj 
fc,m— 1 



^ ^ firnkf^ 



jmk 



k.m—1 



which implies 



F F = nl. 



(19) 



The properties of 9 and + are summarized in the next 
lemma. 

Lemma 1: Let /3, 7 G C s . The mappings 0~ and + 
satisfy 

i. e-(/3) 7 = -e-( 7 )^ 
». e+(/?) 7 = e+( 7 )/3, 

m. 0~(/3)/3 = O, 

i«. e- (e-(/3) 7 ) = [e-(/3),e-( 7 )] ) 

«. e- (e+(/?) 7 ) = {©-(/?), 0+( 7 )}, 

«i. e+ (9-08)7) = [e+C8),e-(7)] = [©-(/?), ©+( 7 )], 

where the commutator of matrices is defined as usual, i.e., 

[A, B] = AB - BA for A,B<E R sxs . 

Proof: Using ( fl3b . one can decompose the left-hand-side of 
@l in terms of the matrices F as 

/ P T F? 7 

0-09)7= : 

V P T F^ 



Every component is the written as 



k,l=l 



y A/iifc7fc 
\ fc,Z=l 



y IkfiklPl 



k,l=l 



y IkfiklPl 
\ k,l=l 



J 



which implies /3 T ^ 7 = --/ T Ff/3. Therefore 8^(^)7 = 
— _ ( 7 )/3. Using dT&b . a similar procedure is applied for 
identity (|m]) in terms of the matrices Dj. The i-th component 
is given by 



/3 T A7 



\ 



k,l=l 



\ k,l=l 



I 



y IkdiuPi 



\ 



k,l=l 



y IkdikiPi 

\ k,l=l 



J 



which gives /? T A7 = 7 T A/3. Thus, 9 + (^) 7 = 6+ (7)/?. 
Identity (|m|) is true since /.yj = for all i and j, and 



y Pifukfik = ^ Pi f uk (3k 

k,l=l 



k,l=l 
k^l 



y PlfilkPk + ^ filfilkPk 
k<l k>l 

^^PlfilkPk — y filfilkPk 



k<l 

=o, 



where the negative sign in the last summand was obtained 
because of the antisymmetry of The left-hand-side of 
identity ([m^ is decomposed as 

e- (e-(/3) 7 ) 




fe=i 



fe=i 



V fc=i 



£ fsskP T F kl 



k=l 



By (lilt , the component of this matrix is 

s 



k=l 

'(- 

3 Tr 



h kFk T 



fc=l 



/3 J [i^-Pfb 
^^7 - /S^-fft 



= p 1 FlFj 1 - 1 1 FiF^ 

= (e-(/3)e-( 7 )-e-( 7 )e-(/3)).. 
= ([e-C8) > e-( 7 )]) ii . 

Similarly, decomposing the (i, j ) component of the left-hand- 
side of ju]) and using (fT~2T > gives 



ftO 



fe=l 



= -P t F 1 D j1 + f3 T DjFij 
= p T F?D n - P T D j FT 1 
= p T F?D jl + 1 T F?D p 

= (e"(/3)e + ( 7 ) + e-(7)e + (/3)) ij 
= ({e-(/3),e+( 7 )}) ir 

Again, decomposing the component of the left-hand- 

side of juT]) and using (TT~3T > gives 



(e+ (e-(^) 7 )) =-Y,d ijk P T F } 



k=l 



— —fi T I E dijkFk 



= -/? T (F l A + FA)7 
= fFjD 3l - -fDiFjp 

= (e-(/?)e+( 7 )-e+( 7 )e-(/3)).. 
= ([e-03),e + ( 7 )]) y . 

Finally, applying the same procedure but using dT4b instead 
gives 



(e+ (e-^)).^-^^/?, 



an 



^=1 



= I l^ijh^h J 7 
= /3 T A'if 7 ~ l T F?D 3 p 

= (e+(/3)e-(7)-e-( 7 )e + (/3)) y . 
= ([e + (/3),e-( 7 )]) u ., 

which completes the proof. ■ 

Lemma 2: Let A, B G C n ™ xs , and denote by A,, i?i their 
respective rows, i = 1, . . . ,n w . Then 

[z, (Ac) T ] = -2i(Q- (Ax )x, ■■■ , Q-{A riw )x) , (20a) 



[x,(Axf}Bx = -2z2(- e_ (^)^ 



fe=i 



e-CAfcje+CBfcjx-ie-CAfcje-^x , (20b) 



{Bx) T [Ax,x T ] = 2iJ^ (-Q-(A k )B 7 k 



fc=i 



+ e-{A k )e + (B k )x + ie-(A k )e-(B k )xj . (20c) 

Proof: The goal is to rewrite ( I20at in terms of [x, x T ] and 
{x,x T } in order to apply ( PT71 ) and ( fl~8b . Then 

[x, (Ar) T ] 

/ xix T Aj ■ ■ ■ xix T A^ 



/ Aixxi 



AiXXi 



\ A nw xxi • • • A Uw xx s 
= (xx T A T 1 - {{A 1 x)x T ) T , • • • , xx T A T Tlm - {{A lX )x T ) T ) 

— — (xx ) j4-^ j * • • j xx -A n ^ i^XX ) i ) 

= ([x,x T L4f,--- ,[x,x T ]<J 

=2i(e-(x)Af,-- - ,8-«). 

Thus, Lemma Q] gives 

[x, (Acf] = -« (e-(Ai)x, • • • , e-(A n „)x) . 



For ( |20bt , note that the scalar operator £>;x commutes with 
Q~(Aj) for any i and j. Recall that 

xx T = l;([x,x T } + {x,x T }). 
It then follows that 

[x, [Axf] Bx 



-2i(e-(Ax)x,--- ,@-(A nw )x) 
AiFfx ■■■ A nw F?x 



= -2i 



A k F?xx T Bl 



2 *E 

fe=1 V A k Fjxx T Bl 

nw ( A k Fl(ll s + e+(x)+iQ-(x))B 



fe=i 



V A k F?(lI s + e+(x)+iQ-(x))Bl 

2i E (-e-(A k )Bl + e-(A k )e+(x)Bl 



fe=l 



■ ie-(^ fc )e-(x)^ 



2i^2 (-Q-(A k )Bl + e-(A k )Q+(B k )x 
■iQ-(A k )Q-(B k )x 



Similarly for d20cb . one has, applying (HJ, that 

(Bx) T [Ax,x T ] 

= 2iJ2 (B k xA k F?x, ■■■ , B k xA k Fjx) 



k=l 



= 2%Y, (B k x 8f {A k )x, ■■■ ,B k x Q- (A k )x) 
k=l 

= 2i£ (Q-(A k )(xx T ) T B T k , • • • , Q-(A k ){xx T ) T B 



k=l 



fc=l 



2i£(e-(A fc )(-/ s + e+(x)-ie-(.T))B, T 



fc=i 



2i£ (-e-(A fc )s, T + e-(A fc )e + (x)i? 



■ <e-(^)e-(x)Bj 



= 2iV [-Q-{A k )Bl + Q-(A k )Q+(B k )x 



fc=l 





and (0) are 




[x,H] = 




[x,£ T ] = 


V B nw x J 


[x,Lt] = 


) ( B nw x ) 










+ 2iO~ 



ie-(A k )e~(B k )x 



The explicit computation of the vector fields in (|2]i is given 
in the next lemma. 

Lemma 3: The component coefficients of equations (O 



(21a) 
(21b) 

(21c) 



k=l 



+ 2iQ-(K*)Q + (K k )x + 2Q-{K*)e-(K k )x\ , (21d) 
-2ie-(A ft )e+(A#)x + 2e-(A fe )e-(A#)xJ . (21e) 



Proof: Commutators (I21ab - (l21cb follow directly from (120a 
Commutator ( I21db is computed out of d20bb as 

[L*,x T ]L = -[x,tf]L 



fe=l 



]T(^e-(Af)A£ + 2ie-(Af)e+(A,) a 



+ 26-(Af)e-(A fc )x^ 
Finally, commutator (121 el l is obtained using d20cb as 

(it [x,L T f) T =-(Lt [i# x r ]) T 



fe=l 



£(^-(a?)a: 



2iQ-(A k )e+(A*)x + 20-(A fc )e-(A#)a; 



IV. Physical Realizability 

The main contribution of the paper are given in this 
section. First, physical realizability is introduced next 

Definition 3: A system described by equations (0 and 
( [Tol l is said to be physically realizable if there exist H and 
L such that © can be written as in (0. 

The explicit form of matrices Aq, A, Bi k , B 2k , C\ and C 2 
in terms of the Hamiltonian and coupling operator is given 
next. 

Theorem 1: Let TL = ax, with a T £ R s , and L = Ax, 
with A e C"™ xs . Then 
4i " 
n 



k=x 

n w 

A = -2Q-(a)+J2(R k -iQ k ), (22b) 



fc=i 



Bik = 0- [i(A* - A fc ) 

s 2fe = e-(A k + A*), 

d=A + A*, 
C 2 = i (A* - A) , 



(22c) 

(22d) 
(22e) 
(22f) 



where 



and 



R k 4 e-(A fe )6-(A#) + e-(A#)6-(A fe ) 

Q fe 4 e-(A fc )e+(A#) - e-(A#)e+(A fc ). 

Proof: The proof follows by direct application of Lemmas Q] 
12 and [3] That is, equation (ffji is re-written using (I21ab - (l21eb 
as the following bilinear QSDE 

dx = -2Q~{a)xdt + —^Q-{A*)Aldt 

11 k=l 

+ g(e-(A fe )e-(Af) + e-(Af)e-(A,))xdt 



k=l 



-iJ2(Q-(A k )Q + (A*)-e-(A*)e + (A k ))xdt 



k=l 



(23) 



+ i(e-(Af-Ai)av ,e-(A* w -A nw )xj dWi 
+ (e~{A 1 +Af)x,--- ,e~{A nw +A*Jx) dW. 

Also, as mentioned in Section [II] the output fields Yi and Y 2 
depend linearly on L, and the input fields W\ and W 2 , 



i.e., 



dYj 
d% 



A + A* 
i(A* - A) 



;dt 



dWi 
dW 2 



It is now easy to identify matrices Aq, A, B\ k , B 2k , Ci and 
C 2 , which ends the proof. ■ 

Note that all matrices involved in the above equation are 
real. To confirm that, observe that A*— A is purely imaginary 
and A + A* is purely real. Now fix k and compute the real 

partof9-(A fe )A£ and 0-(A fe )0+(A# )-0"(A# )6+(A fe ). 
Given that (0~(A fc )Aj.) # = -0-(Afc)A£, one has that 

i?e{6-(A)A t } = 1 (e-(A)At + (0"(A)At)#) = 0. 
Also, 

i? e {9-(A fc )6 + (Af ) - 9-(A# )6+(A fe )} 
= \ (e-(A fe )6+(A#) - e-(A#)0+(A fc ) 

+ (e-(A fc )6+(A#) - 6-(A# )©+(A fe )) # ) = 0. 

Moreover, by direct inspection Bik = —Bi k for 1 = 1,2 
and k = 1, . . . , n w . 

Now, from a control perspective, it is necessary to charac- 
terize when a bilinear QSDE posses underlying Hamiltonian 
and coupling operators which allows to express the matrices 
comprising (0 and dTOb as in Theorem Q] Thus, the second 
and most relevant result of the paper is given in the next the- 
orem, which establishes necessary and sufficient conditions 
for the physical realizability of a bilinear QSDE. 

Theorem 2: System © with output equation dTOb is phys- 
ically realizable if and only if 

i. A = -Y j {iB lk + B 2k ){{C 1 ) k +i(C 2 ) k ) T , 
k=l 

a. Bi fc = e-((c 2 )fc). 
in. B 2k = e-((c 1 ) k ), 

iv. A + A T +J2 B ik B ik T = Je + (A ), 



i,k=l 



where (Ci) k indicates the fc-th row of Cj. In which case, the 
coupling matrix can be identified to be 

A = ^(C\+iC 2 ), 

and a, defining the system Hamiltonian, is 

1 / 1 "™ 

a = -vec [A T A + - £ ([B 2k , + ((C 2 ) fe )] 

V ^ fe=i 



[B lfci e+((d) fc )] F. 



(24) 



Proof: Assuming that (O and ( fTOb are physically realizable 
implies that (I22ab - (l22fb are satisfied. By comparison, condi- 
tions (jllj) - (|m|) hold. Condition (0 is written from (I22eb and 
(l22l as 



A = - V 9-((Ci) fc - x(C 2 ) fe )((Ci) fc +i(C7 2 ) fc ) 3 
71 * — ' 

fe=i 

1 ™» 

= - y)(«Bu + B2fc)((Ci) fc +i(C , 2 ) fc ) r 



fc=l 



Now, one has that 

B lk B lk T = Q-(A* -A k ) 2 

= e-(A#)e-{A*) - e-(A#)e-(A) 
- e-(A)e-(A # ) + e-(A)e-(A). 

Similarly, 

B 2k Bl=-e-(A k + A*) 2 

= e-(A#)e-(A#) - e-(A#)e-(A fc ) 
-e-(A fc )e-(A#)-e-(A fe )e-(A fe ). 

Thus, Bi k Bi k T + B 2k B 2k T = — 2R k . One can now rewrite 
A in terms of a, B lk and B 2k as 

A = - 2Q-(a) 

-^Y,( BikBikT + B2kB2kT )- i ^ k - (25) 

fe=i 

Similarly, 

A T = 20 _ (a) 
i ««» 

-^J2( BikBikT+B2kB2kT )~ i ^ k - (26) 
^ fe=i 

Adding ( f23T > and (f26l > gives 

^4 + A T = - £ (B lk B lk T + B 2k B 2k T ) -i(Q k + Ql) . 

k=l 

The component of Q k + Q\ is computed as 

(Qk + 01)^ = e-(A fe )e+(A#) - e-(Af)e+(A fc ) 
-e+(A#)e-(A fe ) + e+(A fe )e-(Af) 
= - aJf^a* + \:r,n,\, 

-AlD^ + AlD^A*. 

Note that every summand is a scalar, which is equal to its 
transpose. By ( PT3l ) and (TBI) , it follows that 

{Qk + Qk)ij 

= AKDiFj + DjFi) Af + AKFiDj + F j D l )Af 



k=l 



= AlJ2 djikFkA* + AJY, diokFkA 

k=l 

s 

= 2^2d ijk AlD k A* 

k=l 

= 2 (©+(©" (A fe )A#) 



= -2 e+(e-(A*)A fc ) 



= f (e + (^o)) 4J . 



Therefore, adding (1251 1 and ( I26l l gives 



A + A r + £ fl tt B£ = £e+(4>), 



;,fc=i 



which is condition Conversely, one needs to show that 
if conditions (fIj)- (|M>)) of Theorem [2] are satisfied, then there 
exist matrices a and A such that system (0 is physically 



realizable. Let 

e-(a)±\(A T -A 

1 n m \ 

+ 2 E (t^, e + ((C 2 ) fc )] - e + ((d) fc )]) . (27) 
fc=l / 
It is trivial to check that the right-hand-side of ( f27T > is 
antisymmetric and hence this equation uniquely defines a 
via (IT~5T >. Also, let 

A=i(d+iC 2 ). (28) 

Then, Q k can be written in terms of Bi, B 2 , C±, and C 2 as 
follows 

Q fc = i (e-((d +ic 2 ) fe )e+((c 1 - to,)*) 

-e-((d -*c 2 ) fe )0 + ((Ci +«Ca) fc )) 
-e + ((c 1 +ic 2 ) fc )e-((c 1 -«c 2 ) fc ) 
+ e+((c 1 -ic 2 ) fe )e-((Ci+ic 2 ) fc ) 

= -|(e-((Ci) fc )e + ((c 2 ) fc ) 
-e-(((7 2 ) fc )e+((c 1 ) fe )). 

From ([H]) and (|m[). it follows that 

Q k = e-(A fc )e-(Af ) + e-(Af)e-(A fc ) 



= -2(^e+((c 2 ) fe )-i? lfe e + ((C7 1 ) fc )) 



Then. 



Qfc - Q fc = 2 ([^,e+((c 2 ) fc )] - [B lk , e+((Ci) fc )]) . 

Similarly, it is simple to write R k in terms of C\ and C 2 . 
That is, 

i? fe =i (0-((d +*C 2 ) fc )0-((C 1 -<C 2 ) fc ) 

+0-((d -iC 2 ) fe )0-((d +«C a ) fc )) 
= i(0-((C 1 ) fc )0-((C 1 ) fc )-i0-((C 1 ) fe )0-((C 2 ) fc ) 

+ ie-((c 2 ) fc )e-((Ci) fc ) + @-((c 2 ) k )Q-((c 2 ) k ) 
+ Q-((c 1 ) k )e-((c 1 ) k )+iQ-((c 1 ) k )e-((c 2 ) k ) 

-i0-((C 2 ) fe )0-((Cx) fe ) + 0-((C 2 ) fe )0-((C 2 ) fc )) 

= ^ (e-((Ci) fc )e-((cri) fc ) + 0-((c 2 ) fe )©-((c 2 ) fc )) . 

It clear that R k is symmetric. From ([m]) and (|ml) . one obtains 
that 

i? fe = 0-(A fe )0"(A#) + 0-(A#)0"(A fc ) 

= - o { B ^ B Ik + B ^kBj k ) . 



From (0) and 



+ (A o ) = -^(Q fc + or 



Since (JzvJ) implies 



i,fe=l 



one can use ( |27T ) to obtain 

e-(a) 



V fe=i ^ 

+i ([i?2 fc ,e+((C 2 ) fc )] - [B lk ,Q+((C 1 ) k )]) S j\ 



k=l 



1 1 



fc=i 



which is equivalent to ( I22bl i. Moreover, using (fT9] l, (|27] | and 
applying the stacking operator to Q~(a), a is explicitly 
obtained as 



(F T vec(6-(a))) J = (F T Fa T ^ 



not. 



Hence, 



1 / , «» 

a = -vec U T - A + - £[fl 2fc) 9+((C 2 ) fc )] 

V fe=i 



-[B lfc ,e + ((C 1 ) fe )]j F, 

which completes the proof. ■ 

Note that the physical realizability conditions do not require 
the computation of the Hamiltonian (l24b . which depends on 
the structure constants d and /, in order to know whether or 
not the system given by equations (O and (TTOb is quantum. 

V. Conclusions 

A condition for physical realizability was given for open 
multi-level quantum systems. Under this condition it was 



shown that there exist operators TL and L such that the 
bilinear QSDE (O with output equation ([Tol l can be written 
as in (f2]i. This condition used explicitly the algebra generated 
by SU(n). Moreover, the interaction of the system with 
multiple quantum fields was introduced to the formalism. 
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